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Abstract 

In almost all structure formation models, primordial perturbations are created within a homoge- 
neous and isotropic universe, like the one we observe. Because their ensemble averages inherit the 
symmetries of the spacetime in which they are seeded, cosmological perturbations then happen to 
be statistically isotropic and homogeneous. Certain anomalies in the cosmic microwave background 
on the other hand suggest that perturbations do not satisfy these statistical properties, thereby 
challenging perhaps our understanding of structure formation. In this article we relax this tension. 
We show that if the universe contains an appropriate triad of scalar fields with spatially constant 
but non-zero gradients, it is possible to generate statistically anisotropic and inhomogeneous pri- 
mordial perturbations, even though the energy momentum tensor of the triad itself is invariant 
under translations and rotations. 
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I. INTRODUCTION 



The cosmological principle states that the universe is homogeneous and isotropic on large 
scales. Although it was originally invoked for the sake of simplicity, it is by now solidly 
grounded by observations. A Friedman-Robertson- Walker metric accurately describes the 
evolution of the universe on large scales, and it suffices to explain phenomena that range 
from the abundance of the light elements to the (near) isotropy of the cosmic microwave 
background. However, homogeneity and isotropy are clearly just approximations. Because 
of the same cosmological principle, we are able to describe inhomogeneities and anisotropies 
in our universe as small perturbations around the lowest order homogeneous and isotropic 
background. 

There is a generalized sense of homogeneity and isotropy that may also apply to these 
perturbations, namely, statistical isotropy and homogeneity. We shall not extensively dwell 
on these concepts, and simply note that perturbations are statistically homogeneous and 
isotropic if their ensemble averages are invariant under translations and rotations (see for 
instance pp.) Given that the cosmological principle appears to be solidly grounded by obser- 
vations at this point, it is time to inquire — instead of postulating — whether perturbations 
also satisfy the generalized, statistical, cosmological principle. Certain anomalies in the tem- 
perature anisotropies in the CMB actually suggest that the answer is negative [H [5l [6l [7] . 

If perturbations and the underlying background were completely independent from each 
other, the symmetries of the latter would not have to bear any relation with the statistical 
symmetries of the former. However, this independence is not in agreement with our ideas 
about the origin of structure. According to our present understanding, the structures that we 
observe originated from tiny primordial perturbations that were generated during a seeding 
stage in the early universe [21 E]. As a result, primordial perturbations always seem to 
carry the symmetries of the background on top of which they were seeded. Inflation for 
instance is typically driven by a homogeneous scalar field, and as consequence, primordial 
perturbations created during infiation are statistically homogeneous and isotropic. It is 
also possible to seed primordial perturbations in an anisotropic but homogeneous universe, 
and, in agreement with our claim, the resulting primordial perturbations turn out to be 
statistically homogeneous but anisotropic [HI [9]. If the universe is inhomogeneous during 
the seeding stage, the statistical properties of primordial perturbations will accordingly 
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reflect those spatial variations [TO] . 

The conclusion that seems to emerge from our previous discussion is that if the universe 
is homogeneous and isotropic during the seeding state, primordial perturbations have to be 
statistically homogeneous and isotropic. In this article we argue that this conclusion does not 
necessarily hold. We present a class of models for the origin of structure in which the universe 
is exactly homogeneous and isotropic throughout the seeding stage, with the peculiarity that 
the resulting primordial perturbations are statistically anisotropic or inhomogeneous. 



II. A SCALAR TRIAD 



The class of models that we are going to consider bears a significant relation to previous 
works in the literature. It contains a set of three arbitrary non-canonical scalar fields ^Jq, 
(where a = 1,2, 3) minimally coupled to Einstein gravity pT| [T^. The action of this "scalar 
triad" is 



S = d Xa/— ^ 



where the Lq, are functions of a single argument^ 



(1) 



Xa, = ]^g^''d^'^c,d^'^^, (2) 

and Cm is the Lagrangian of the remaining matter fields ip, which can include couplings to 
the triad. From a different perspective, the triad can be regarded as part of the gravitational 
sector; the action ([T]) then leads to what one could interpret as a modification of gravity 
|13j . Note that throughout this work we use Einstein's summation convention; a sum over 
indices is implied only if they are in opposite positions. Otherwise, the location of internal 
and spatial indices does not have any significance. 

The stability of the system described by the action ([T| imposes conditions on the functions 
La- The requirement that perturbations around any background value of have a positive 
definite Hamiltonian (in order for the theory to be ghost-free) implies [H] 



^ To simplify comparisons with the Hterature we work with a (+, — , — , — ) metric signature. 
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In turn, the last two conditions guarantee that the system has a well-defined initial value 
formulation, at least locally, that is, that scalar perturbations propagate around appropriate 
characteristic cones. 

The field equations of motion take the form of a conservation law, which is an expression 
of the symmetry of the action under shifts of the scalars, 

dJv^^9'-'d.^e]=0. (4) 



The energy-momentum tensor of the triad is simply the sum of the energy-momentum tensors 
of the individual scalars, 

III. BACKGROUND 

We are interested in cosmological solutions that satisfy the cosmological principle, that 
is, in spatially homogeneous and isotropic universes, 

ds'^ = a^{ri) ■ (dri^ - dx'^). (6) 

As opposed to what one would naively expect, homogeneous configurations do not require 
the field itself to be homogeneous. Since the action ([T]) only depends on the three scalars 
through their derivatives, the actual value of the fields is physically irrelevant, since only 
their gradients have physical meaning. Therefore, homogeneous configurations only require 
that the field gradients be spatially constant. Hence, rather than assuming that the spatial 
gradients are zero, as is typically done, we shall study the ansatz 

</^a = A'x", (7) 

which in conjunction with equation (|6| results in 



1 A"^ 

A-.X„ = ---. (8) 



Note that we are implicitly assuming that A is a constant. Indeed, one can readily verify 
that for time-independent A, the ansatz ([T]) solves the equation of motion (|4]). As a matter 
of fact, for the metric ^ the equation of motion reads 



dXa 



(^^'L + a^y^; - Ay^.) + ^ ■ - dx^ ■ d^c) = o, (9) 



which is clearly satisfied by the ansatz ([T]). Note that a prime denotes a derivative with 
respect to conformal time, and A = did^ is the fiat space Laplacian. 

Although the configuration ([T]) does respect homogeneity, it does appear to violate 
isotropy, since the field gradients diipa = 5ia are not invariant under spatial rotations. 
However, this impression is misleading, because the contributions of each scalar to the en- 
ergy momentum tensor (which is all gravity sees) may add to an isotropic energy momentum 
tensor. In fact, the total energy momentum tensor of the three scalars is given by 



T\ = (10a) 

a 

T\ = 0, (10b) 



T'j = 2Xa---^5\6ja - LgS'j. (10c) 



Therefore, if dLa/ dlogX^ does not depend on a, the energy-momentum tensor is invariant 
under spatial rotations. The simplest way to satisfy this condition is to assume that the 
Lagrangians are identical, 

L = L^ = Lp {a^(3). (11) 

In this case the theory has a Z3 symmetry that acts on the three scalar fields by permutation. 
For Lagrangians linear in X, this Z3 is just a subgroup of 0(3). 



If conditions (11) are satisfied, the stresses are isotropic, and we can therefore identify 



P = -3L, (12) 
dL 

p = 3L-2X- (13) 

as, respectively, the energy density and pressure of the scalar triad. In this sense, the three 
scalar fields are analogous to the group of three vectors described in [15] ; the field gradients 
play then the role of of the spacelike vectors. Let us note that the ghost-free condition 
dL/dX > implies, in conjunction with X < 0, that p + p > 0. The scalar triad cannot 
behave like a phantom [TB]. 

For the sake of illustration, it is going to be convenient to have a class of concrete 
Lagrangians in mind. For our purposes it will suffice to consider a simple but sufficiently 
general set of models with Lagrangian 

L = C- {DXY, (14) 
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where C and D are two arbitrary dimensionful constants and n is an arbitrary exponent. In 
this case, the triad always behaves hke a fluid with equation of state 

2ri-3 

w = (15) 
The Lagrangian satisfies the stability conditions ([s]) if n > 1/2. 



IV. STATISTICAL ANISOTROPIES 



Our first goal is to show that within the class of models that we are studying it is possible 
to create statistically anisotropic scalar primordial perturbations. We shall first define what 
we mean by "scalar" and show that in our decomposition of the perturbations scalars, vectors 
and tensors generally couple to each other. 



A. Perturbation Decomposition 

Let us start by considering metric perturbations in longitudinal gauge, 

(is^ = [(1 + 2(^)dr]'^ + 2Si dr] dx' - {6,j - 2^ % - hij)dx'dx^] , (16) 

where Si is a transverse vector, and hij is a symmetric transverse and traceless tensor. 
Essentially, we just have decomposed the metric perturbations in irreducible representations 
of a subgroup of symmetries of the FRW metric, namely, spatial rotations. Rather than 
proceeding similarly with the triad perturbations, it will be instructive to expand^ 

S^o.=A^-{d^s + Vo,), (17) 

where s is a scalar and Va is a transverse vector, daV"" = J2a ^a^a = 0. Indeed, our scalar 
field background ([T]) is not invariant under spatial rotations, but instead, remains unchanged 
under a 0(3) group of symmetries 

Xi^R^jX^ ^o.^ R-'pifi", RE 0(3). (18) 



The decomposition (17), as well as (16), is a decomposition of perturbations in irreducible 



representations of the latter group of transformations. 



^ I thank Ignacio Navarro for pointing this out. 
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We shall not write down the different components of the linearized Einstein tensor, which 
can be found in many different monographs on the subject (e.g. [T7]) . The components of 
the linearized energy momentum tensor are 



6T\ 



6T' 



aX 

-2X^{S, + d,s' + v',) 
dL 



6Tj = 2X— ■ [^6,j + {2did, - 5,jA)s + d^vj + d^v, + Kj] + 
+ 2X2^ ^ 5a.5a, [2^- + 2dls + 2d^v^ + . 



(19a) 
(19b) 

(19c) 



Although equations (19a) and (19b) appear to suggest that scalar, vector and tensor pertur- 



bations decouple, closer inspection of equation (19c) shows that this is not the case whenever 
(PL/dX"^ 7^ 0. Consider for instance the decomposition in irreducible representations of the 



last term in equation (19c) 



5ai5aj haa = A5ij + didjB + diCj + djCi + Dij. 



(20) 



The scalar A is given hj A = — ^ J2a ^aa [E], which does not vanish and does depend on 
the tensor hij. Hence, the linearized Einstein equations couple scalars to vectors and tensors. 
The same conclusion holds if the triad perturbations are expressed in terms of the variables 
S(poi- The violation of the decomposition theorem was noted in a similar context in 



Note by the way that the term 2XdL/dXhij in equation (19c) gives a mass to the graviton 

m. 



There are two cases where scalars, vectors and tensors can be decoupled, namely, when 
X'^d^L/dX'^ = or when the metric perturbations themselves are negligible. In the first 



case, the group of transformations (18) is not just a symmetry of the background, but also 



a symmetry of the action. As a result of this enhanced symmetry the triad perturbations 
turn out to be statistically isotropic, as we shall see below. We are hence led to the limit in 
which the metric perturbations are much bigger than the triad perturbations, and are thus 
neghgible. We expect this to be the case, for example, if the expansion of the universe is 
sufficiently close to de Sitter. 
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B. Primordial Perturbations 



We shall study now the perturbations in the triad created during a stage of seeding in the 
early universe. At this point it is not important to specify what is the dominant component 
of the universe during that stage; it could be the triad itself, or it could be a different scalar, 
like in curvaton models [20]. We shall simply assume that the early universe inflates, and, 
for definiteness, that its effective equation of state w is constant, 

/ r/ V 2 
a = oo — , where p = . (21) 

\VoJ l + 3w 

In order to calculate the primordial spectrum of metric perturbations, we have to specify 
some of the details of the transition between the seeding stage and a radiation dominated 
universe. If reheating is instantaneous and occurs along an hypersurface of constant Xiamen 
we find that after a stage of inflation with equation of state w the universe contains metric 
perturbations (in longitudinal gauge) even though they were negligible during inflation |2T] , 

3w — 1 ^ , ^ \ ^ 1 6Xa 1 \ ^ „ ^ 
<^^—^S, where ^ = }_^-— = j^}^dj^^. (22) 

a a 

Similarly, if the decay rate of the component that drives the seeding stage depends on the 
squared gradients of the triad Xa, perturbations in the latter will be converted into adiabatic 
primordial perturbations at the end of the seeding stage [22] . In our theory, by assumption, 
the Lagrangian is symmetric under shifts of the scalars and permutations, so the decay 
rate has to be of the form Ftot = J2a^(-^a)- As a consequence, we expect primordial 
perturbations of the order 

^^i5iV.^^diogr 

9 Ttot 27rflogX ' ^ ^ 



where 6 has been defined in equation (22). The term proportional to log F/rf log X is 
typically of order one, an hence, primordial perturbations are again proportional to 6. In 
the following we shall calculate the spectrum of S. 



C. Power Spectra 

We are ready now to calculate the action for the cosmological perturbations. As men- 
tioned above, we are neglecting metric perturbations, so all we have to do is expand the 
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action of the scalars to quadratic order in the presence of an unperturbed FRW metric. 
Proceeding in this manner we arrive at 



5S = ^ j (fx 



dL 
dX 



(24) 



from which one can directly derive the stabihty conditions (|3]). Note that we have expressed 
the field perturbations in terms of the perturbed fields 6ipa- These are the fields that decouple 



in the absence of metric perturbations. When we substitute equation (17) into the action 
(24), we find that the term proportional to d'^L/dX'^ couples the vector Va to the scalar 



s. The very same term manifestly breaks rotational invariance in the action, and will be 
ultimately responsible for the generation of statistically anisotropic perturbations. 

In order to quantize the perturbations, it is convenient to further decompose them in 
Fourier modes of an appropriately rescaled scalar field (with periodic boundary conditions), 

o dL 



b ■ 5ipc 



where h 



dX' 



(25) 



Note that h plays the role of an "effective" scale factor, sometimes referred to as the pump 



field. Substituting the expansion (25) into the action (24), and employing the particular 



class of models in ( 14 ) we arrive at 



5S^ = VY, j dv\ 
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h" 



VkV-k 



where 



k'' + 2{n-l)kl 



(26) 



(27) 



is the anisotropic dispersion relation of the field perturbation 5ipa and V is the comoving 
volume of the universe.. As the reader may have expected, the dispersion relation is isotropic 
only if d'^L/dX'^ vanishes (ra = 1). 



The quantization of (|26|) proceeds along the conventional lines [23]. The properly nor- 

(28) 



malized mode functions are 



— TlTj 



where is the Bessel function of the first kind and 

1 - 2(2 - n)p 



(29) 



Using the definition of 5 in equation (22 ), and the relation between Va and Sipa in equation 



(25) it is now straightforward to calculate the spectrum of 6. For z/ > the spectrum is 



constant on large scales {ujaf] <^ 1). Therefore, in that limit it can be evaluated at any 
convenient time, which we choose to be the time at which a reference mode k^, crosses the 
horizon, k^/a^ = H^. On large scales the primordial spectrum is then 



where the dispersion relation uj^ is given by equation (27), and the index v is given by 



equation (29). The equation of state of the dominant component during inflation is w, and 
p* and denote, respectively, the energy density and pressure of the triad at the time of 
crossing. 

If n 7^ 1, primordial perturbations are statistically anisotropic, even though the back- 
ground energy-momentum tensor and the spacetime metric are invariant under rotations. 
In any case, the primordial spectrum still possesses certain remnant symmetries. In partic- 
ular, it is scale free and invariant under spatial inversions, k —>■ —k. The scaling symmetry 
follows from the structure of the power spectrum, which is a product of two functions, 
V{k) = V{k) ■ VL{k), where k = k/k. One could ask whether this factorization is generic, 
or whether it is just particular feature of our set of models. By tracing back the steps we 
took to derive the power spectrum, we see that the origin of this property is the form of the 
dispersion relation Ua- If uja consists of terms with different powers of ka, Vsik) cannot be 
factorized. Similarly, it is easy to see that the parity symmetry arises from the symmetry 
of the dispersion relation. Dispersion relations with odd powers of k would lead to parity 
violations. 

As we shall argue below, the factorization of the primordial spectrum allows us to identify 
an effective spectral index of the perturbations, — 1 = 5 — 2z/. In the following we shall 
focus on the phenomeno logically preferred case of scale invariant perturbations, = 1. The 
spectrum is scale invariant in this generalized sense if u = 5/2, which implies 

n — 3 , , 

w = — . (31) 

We can compare this condition with the equation that relates the expansion of the universe 
to the parameter n when the triad dominates the energy density of the universe, equation 



(15). Only for n = are both mutually consistent. However, for this value of the parameter 
perturbations are quantum-mechanically unstable. Hence, in this class of models, the com- 
ponent that drives inflation cannot be the triad. Because stability requires n > 1/2, and 
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since in our context the seeding requires a stage of inflation, —1 < w < —1/3, the range 
of values of exponents compatible with a scale invariant spectrum is 1/2 < n < 2. These 



conclusions are closely tied to our class of Lagrangians (14), and could be relaxed in other 
types of models. 

D. Observational Signatures 



Let us briefly describe what type of signatures a spectrum of the form (30) would leave 
on the cosmic microwave background. The impact of an arbitrary statistically anisotropic 
spectrum of primordial perturbations on the CMB was calculated in pTj, so all we have to 
do is apply the formalism developed in the last reference, which the reader may consult for 
further details. 

Our starting point is the decomposition of the primordial spectrum in spherical harmon- 
ics. 

The functions Vem quantify the power in the harmonic im as a function of scale k. For 
statistically isotropic perturbations, only the scalar component "Pqo differs from zero. Cer- 



tainly, this is not the case in our class of models. In particular, comparing equation (30) 



with equation (32), and remembering that the spectra of S and $ just differ by an overall 



constant, we arrive at 



'P^rn = ^- [ d'kY.Uk) . J2 7 V^' (33) 



where A is an overall amplitude we shall not worry about. Because all the harmonics share 
the same power of k we can identify Hs = 6 — 2u as the spectral index of the primordial 
spectrum. Since the spectrum is invariant under spatial inversions, harmonics with an odd 
angular momentum quantum number vanish. Vim = for odd i. Similarly, because Q[k) is 
invariant under reflections on the k^ — ky plane, Vem = for odd m. 

The different harmonics of the power spectrum are connected to the statistical properties 
of the CMB. In particular, the expectation value of the estimator typically used to measure 
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FIG. 1: A density plot of the function V^{k) for n = 3/2. 



the angular power spectrum^ is 

Ce = ^ = / f ^'(^) ^oo(^)' (34) 

m 

where is the conventional transfer function ^25j . Therefore, the angular power spectrum, 
as defined above, does not contain any information about statistical anisotropics, since it 
is only sensitive to the scalar component of the power spectrum Poo- If the spectrum is 
scale invariant (in our generalized sense), statistically anisotropic perturbations will fit the 
angular power spectrum as well as statistically isotropic ones. 

The anisotropic components of the power spectrum affect the non-diagonal two-point 
functions, 

= (-^)'^^'^ $^/^(^i,mi;^,m|4,m2) j yA,,(fc)A,,(A;)P,^(fc), (35) 

where Z) is a product of Clebsch-Gordan coefficients [1]. Hence, statistically anisotropic 
perturbations yield correlations between the different spherical coefficients of the CMB. A 
color-coded contour plot of the power spectrum P$ on the unit sphere {k G S'^) is shown in 

In practice, sky cuts and noise in the temperature maps require the use of somewhat different estimators 
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TABLE I: The first no n- vanishing multipole components of the power spectrum for n = 3/2. 



figure [l| From the figure, it is apparent that the smallest structures in P<i>(fc) have angular 
sizes of about tt/6. Hence, we expect its multipole components to be insignificant beyond 
i ^ 8. In table [T] we also list the first few non-vanishing multipoles of our fiducial model 
with n = 3/2. Note that there is no dipole, which is the multipole we would expect from 
just one preferred spatial direction. Instead, the dominant multipole is £ = 4. 



Statistically anisotropic primordial perturbations of this type could a priori explain the 
anomalous alignment of the quadrupole and the octopole of the CMB [4J. Although no 
evidence for statistical anisotropy in the CMB has been found in the model-independent 
analysis of , this does not imply that statistically anisotropic primordial perturbations are 
ruled out. Bounds on the amplitude of the different power spectrum components Vim do not 
exist (as far as the author knows), so one cannot exclude this type of perturbations without 
further analysis. The situation here resembles to some extent the status of non-adiabatic 
perturbations. Although there is no evidence for non-adiabatic primordial perturbations, 
current bounds still allow significant isocurvature components p7] . 

V. STATISTICAL INHOMOGENEITY 

In the light of these results, the reader may wonder whether it is also possible to seed 
statistically inhomogeneous primordial perturbations in a Friedman-Robertson- Walker uni- 
verse. Indeed, it is possible to do so by dropping the requirement that the triad be invariant 
under shifts. Consider for instance the following matter Lagrangian, 

1 1 ^ 

Cm = ^d^xd^X- ^^m^{Va)x^ + L.^{g^^,i),x)- (36) 
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It contains an additional scalar x, whose field-dependent mass breaks the shift symmetry 
of the theory. The Z3 symmetry under permutations of the triad is still preserved, and this 
time we assume that only the scalar field x has couplings to the standard model matter 
fields if). Because the shift symmetry is broken, radiative corrections should lead to many 
additional field-dependent terms in the Lagrangian. Their absence is likely to require certain 
degree of fine-tuning. 

We shall assume that the field x is at the bottom of its potential, x = 0- This is not only 
consistent, but also expected if the mass of the field is sufficiently big. Then, the scalar x 
does not contribute to the energy density of the universe, and does not affect the equations 
of motion of the triad. In particular, the configuration ([t]) is still a cosmological solution, 
and equation (Is]) still gives its energy-momentum. 



A. Power Spectra 



Let us study now what type of primordial perturbations could be produced under such 
circumstances. First, note that the field x couples to matter, so perturbations in the former 
can be converted into adiabatic primordial perturbations at the end of a seeding stage by 



the Dvali-Gruzinov-Zaldarriaga mechanism p2J we discussed in Section IV B Hence, it will 
suffice to calculate the power spectrum of x- As before, we shall neglect metric perturbations. 

In order to obtain the action for the perturbations of Xi l^t us work again in the Fourier 
space of the rescaled variable v = ax- In terms of this variable, the action for the scalar x is 



kX 



V 
2" 



kl 



v^,v_i - k Vk v^i H Vk v-i] +a mi__^ Vk v^i 

a 



(37) 



where ml are the Fourier modes of the squared mass of x evaluated along the configuration 
of the triad fields, 



nil. 



1 

V 



(38) 



which do not depend on time. The form of the action manifestly shows that the Fourier 
components of the field x do not decouple. Consider for instance the equation of motion of 
the mode fc. 



v'l + \ k 



Vk + ^a^ ml_[ Vi = 0. 



(39) 
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It only contains that mode if m1_i oc Ski, that is, if is spatially constant. But in that case 
the action is invariant under spatial translations, so it is not possible to seed statistically 
inhomogeneous perturbations. We therefore encounter again a link between violations of 
the decomposition principle and the generation of statistically asymmetric perturbations. 

Because different Fourier modes are not decoupled, the quantization of this theory requires 
more effort. Instead of diagonalizing the action, we shall treat the mass as a perturbation. 
This is a good approximation for a sufficiently light field, m^/H'^ ^ 1, and a sufficiently 
small number of e-folds of inflation [2H]- To first order in perturbation theory, a non-zero 
mass term generates a correction to the zeroth order power spectrum 



where ifmt is the interaction Hamiltonian, 

^irit = \^<^^ ^tk Vk V-l. (41) 

kl 

Note that we are working in the interaction picture, where the operators Vk carry the free 
(massless) time evolution. 

For simplicity, we shall assume a de Sitter expansion of the universe, p = —1, in which 
the properly normalized mode functions are 



Then, at late times, the integral in equation (40) diverges with the logarithm of the scale 
factor [28]. We shall restrict our calculation to the leading logarithmic correction to the 
correlation function in the limit of late times. Isolating that contribution, after a fair amount 
of calculation with exponential integrals, we find that the two-point function of x is given 
by 

V-k\ . , H\ P /I 1\ 2 , k + l 

■(xlxi) ^ -r-,5ki + ■ + « Wtfc log-TF- (43) 



Hence, the two-point function is not proportional to a momentum-conserving Kronecker 
delta, a momentum dependence that signals the breakdown of statistical homogeneity [1]. 
In this case, the conventional notion of "power spectrum" V^{k) is not defined. 



15 



VI. SUMMARY AND CONCLUSIONS 



In this article we have explored the relation between the symmetries of a Friedman- 
Robertson- Walker universe and the statistical symmetries of its perturbations. Our starting 
point has been the cosmological principle, the homogeneity and isotropy of the universe, 
which we have taken for granted. We have then explicitly shown that it is possible to 
generate statistically anisotropic and inhomogcncous primordial perturbations if the early 
universe contains a scalar triad. The triad is a set of three scalar fields with non-vanishing 
spatial gradients. Its properties are essentially determined by the cosmological symmetries, 
which require that the gradients be mutually orthogonal and spatially constant. Although 
this configuration is not isotropic, it turns out that its energy momentum tensor is invariant 
under rotations, and thus, is compatible with the symmetries of a homogeneous and isotropic 
FRW universe. 

In our seeding scenario, statistically anisotropic or inhomogeneous primordial perturba- 
tions arise from the breaking of rotational or translational invariance by the triad configura- 
tion, which also prevents the decoupling of the different modes of cosmological perturbations. 
Because rotational and translational invariance must be ultimately broken in any model that 
singles out particular directions or locations, it is likely that the connection between sta- 
tistically anisotropic or inhomogeneous perturbations on one hand and violations of the 
decomposition theorem on the other is generic. 

We have discussed a particular class of simple models that lead to a statistically 
anisotropic primordial spectrum of perturbations. Even though in general it is not possible 
to do so, in our example we can define a single spectral index that uniquely captures the scale 
dependence of the primordial power. Requiring scale invariance of the created perturbations 
then imposes additional constraints on the scalar field Lagrangian, and implies — at least for 
the restricted class of Lagrangians that we have studied — that the scalar triad cannot be the 
dominant component during the seeding stage in the early universe. The resulting angular 
dependence of the primordial spectrum is quite rich, and does not simply correspond to the 
dipole pattern one would associate to a single preferred direction. 

It is also possible to generate statistically inhomogeneous perturbations during a stage 
of infiation. In this context, the breaking of translation invariance arises from the spatial 
variations of the triad, and the explicit breaking of a shift symmetry by the coupling of 
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the triad to matter. Because the definition of power spectrum imphcitly assumes statistical 
homogeneity, we cannot quote a power spectrum in the conventional sense, since the two- 
point function now depends on two momentum variables. In the class of models we have 
considered, the kinetic terms of the triad did not play any role, and could be taken to be 
canonical. 

Whether or not primordial perturbations are statistically homogeneous and isotropic is 
a question that should be settled by experiment, not by theoretical prejudice. Certain 
anomalies of the CMB temperature anisotropics indeed suggest for instance that cosmo- 
logical perturbations are statistically anisotropic. Although these anomalies could have a 
primordial origin, we have not explored this possibility here. Instead, we have simply argued 
that statistically anisotropic or inhomogeneous perturbations do not necessarily clash with 
our understanding of the origin of structure in the universe. 
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